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1.  INTRODUCTION. 

Recent  advances  in  microcomputer  technology  have  generated  renewed 

interest  in  the  design  and  implementation  of  sampled-data  control  systems, 

1  2 

a  subject  area  with  an  already  mature  and  well -developed  theory.  *  The 
introduction  of  digital  flight  control  systems  into  modern  aircraft,  for 
example,,  has  for  the  past  several  years  been  receiving  considerable  atten- 
tion.  ”  The  compactness,  flexibility,  accuracy,  and  reliability  of 
flight  control  computers  make  them  ideally  suited  for  the  exacting  guidance 
and  control  requirements  of  high  speed  flight  environments.  Desired  future 
full-scale  integration  of  the  diverse  control  functions  of  aircraft  sub¬ 
systems  will  necessitate  efficient  utilization  of  digital  processors. 

Accordingly,  this  discussion  focuses  upon  a  simple,  direct,  and  computa¬ 
tionally  efficient  signal  tracking  approach  to  the  digital  control  of 
multivariable  linear  systems.  ^ 

Signal  tracking  has,  of  course,  provided  an  effective  design  criterion 
for  a  variety  of  practical  control  problems.  The  popular  linear  quadractic 
regulator^  permits  the  close  tracking  of  reference  variables  (set  points) 
which  are  constant  over  relatively  long  periods  of  time  while  at  the  same 
time  constraining  the  amplitudes  of  control  inputs.  The  MRA  techniques 

g 

surveyed  by  Landau  have  as  their  common  objective  the  minimization  of 

deviations  between  reference  and  controlled  system  outputs  through  the 

adjustment  of  system  parameters  and  various  gain  matrices.  Kuo,  Singh, 

9  10 

and  Yackel  ’  have  employed  output  matching  at  sampling  instants  to  digita- 
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lize  the  control  of  continuous  feedback  systems.  Mehra  *  ’  has  used 

output  matching,  in  conjunction  with  an  impulse  response  plant  representation, 

to  drive  system  outputs  to  desired  set  points  along  constructed  nonlinear 
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This  paper  presents  a  state  space  output  matching  approach  to  the 
digital  control  of  continuous  multivariable  linear  systems.  The  control 
strategy  adopted  relies  simply  upon  the  sequential  generation  of  digital 
inputs  which  force  system  outputs  to  closely  track  prescribed  trajectories. 
The  problem  formulation  is  quite  general  and  is  very  similar  in  flavor  to 
that  of  Mehra.  It  provides,  as  will  be  seen,  an  alternate  approach  to  the 
problems  considered  in  references  4  and  9.  The  control  technique  developed 
relies  upon  only  the  most  elementary  principles  of  linear  systems  theory 
and  hence  is  characterized  by  its  conceptual  and  computational  simplicity. 


2.  THE  CONTROL  PROBLEM. 

Consider  the  sampled-data  feedback  control,  as  depicted  in  Figure  1, 
of  the  time- invariant  dynamical  system 


i(t)  3  Ax(t)  ♦  Bu(t) 

y(t)  3  Cx(t)  +  Du(t) ,  x(0)  =  xQ. 

Here  T  >  0  is  the  sampling  period,  x(t)  e  Rn,  y(t)  eRl,  z(t)  fill*  , 


(1) 


Figure  l  -  Sampled-data  Control  System 


r(t)  e  R*  ,  and  u(t)  e  R  satisfies 

u(t)  =  u(kT),  kT  <  t  < (k+l)T,  k  =  0,1,2...  (2] 

The  matrices  A,  B,  C,  and  D  are  of  appropriate  dimensions.  For  later  con¬ 
venience,  if  all  states  are  observable,  assume  that  the  observer  is  omitted 
and  that  x(kT)  is  fed  back  to  the  controller. 

The  function  z(t)  represents  a  desired  trajectory  that  the  plant  output 
y(t)  is  to  track  closely.  It  may  be  the  output  of  a  dynamical  system  as 
predicted  by  its  model  driven  by  r(t),  or  some  other  appropriately  generated 
signal.  However  z(t)  is  specified,  it  is  assumed  that  z((k+l)T)  is  known 
at  t  *  kT  (See  Figure  1).  Close  tracking  is  insured  here,  as  in  references 
4,  S,  9-11,  by  determining  a  control  variable  sequence  u(0),  u(T),  u(2T),... 
that  insures  close  matching  of  the  components  y^(t)  and  z^(t)  at  t  *  kT, 
k  »  0,1,2,... 

As  will  be  seen,  in  the  absence  of  control  constraints,  the  solution 
of  the  control  problem  considered  reduces  simply  to  determing  optimal  con¬ 
stant  gain  matrices  F  and  G  and  a  sequence  of  vectors  c(kT)  for  the  control 
system  of  Figure  2.  This,  of  course,  is  the  attractive  form  of  a  standard 
linear  feedback  control  system. 


Figure  2  -  Constant  Gain  System 
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3.  OPTIMAL  CONTROL  STRATEGIES. 

Since  controls  u(t)  in  (2)  are  piecewise  constant,  solutions  x(t)  of 
(1)  satisfy  the  difference  equations 


x((k+l)T)  =  $x(kT)  ♦  Qu(kT),  k  *  0,1,2,...  , 


(3) 


where 


$  =  e‘ 


AT 


and  0  = 


f 

'0 


eACT'^  B  dr 


(4) 


Using  these  equations,  a  precise  formulation  of  the  output  matching  control 
problem  of  section  2  is  easily  given.  Specifically,  at  each  sampling 
instant  t  *  kT,  k  >  0,  given  y(kT)  and  z[(k+l)T),  consider  the  unconstrained 
minimization  of 


JCu(kT)) 


*  4  l  q-tz.CCk+DT)  -  y^Ck+UT)]2 

i»l 

1  £  2 
■  2  l  q,[t,((k>l)T)  -  C.$x(kT)  -  (C.O+D. )  u(kT)  ] Z  . 

i*l  1  1  1  11 


CS) 


Here  q^  >  0,  l  <  i  <  l,  Ci  and  are  the  ith  rows  of  C  and  D  respectively, 

and  x((k+l)T)  is  given  by  (3).  The  control  u(t)  determined  through  (2)  by 

the  variables  u(0},  u(T),  u(2T),...  which  minimize  the  functions  J(u(kT)) 

sequentially,  k  *  0,1,2,...,  will  be  said  to  be  optimal.  If  the  more  typical 

problem  of  minimizing  £  J(u(kT))  is  considered,  then  the  strategy  considered 

k>0 

here  clearly  may  be  suboptimal.  However,  if  the  plant  is  sufficiently 
responsive  and  T  is  sufficiently  small,  the  minimization  of  J(u(kT))  should 
provide  desired  close  tracking  of  z(t). 


s 


Consider  the  minimization  of  J(u(kT))  when  all  states  are  observable 
and  controls  are  unconstrained  (Figure  1  with  the  observer  omitted  and  x(kT) 
fed  back).  In  theory,  perfect  output  matching  (J(u(kT))  *  0  for  all  k) 
can  be  achieved  if  and  only  if  the  equations 

(C0  +  D)u(kT)  *  z((k+l)T)  -  C$x(kT)  (6) 

g 

have  a  solution  (cf.  Kuo,  Singh,  and  Yackel  ).  If  a  solution  exists  and  if 
rank  (C9  ♦  D)  *  r,  then  there  exists  a  permutation  a  of  the  indices 
U,2,...,m},with  a(i)  <  a(i+l),  and  matrices  M  and  N  such  that 


rwkT>i 

l*J  -  M[z((k+1)T)  -  c*x(kT)-N 

ua(r+l)(kT) 

LW)CkT)-l 

-uaoio(kT)  - 

This  follows  from  the  row  reduction  of  r  linearly  independent  rows  of 
G0  +  D.  For  r  <  m,  N  consists  of  columns  a(r*l),  a(r+2) , . . . ,  a(m)  of  C0  +  D. 
The  values  of  the  variables  ua(i)(kT),  r  ♦  1  <  i  <  m,  may  be  conveniently 
assigned  arbitrarily.  Thus,  the  solution  vector  u(kT)  may  be  written 

u(kT)  *  F  z((k+l)T)  -  Gx(kT)  -  c(kT).  (8) 

The  rows  F0^j  of  F  are  given  by  for  1  <  i  <  r,  and  Fa^  *  0 

otherwise.  The  matrix  G  is  given  by  G  »  FC$  and,  letting 

T 

c' (kT)  *  FN[u0£r+1j(kT),... »ua^m)(kT)]  ,  the  correction  vector  c(kT)  has 

ca(i)(kTi  *  ca(i)(kT)*  1  f  1  5  r»  and  ca(i)(kT)  3  'ua(i)(kT5  otherwise 
(superscript  T  denotes  transpose).  If  r  ■  m,  F*K*  (C0  ♦  D)"1  and  c(kT)  2  0. 

Importantly,  F  and  G  (and  effectively  the  vectors  c(kT),  since  ua(ij(kT)» 
r  ♦  l  <  i  <  m,  is  arbitrary)  are  independent  of  k.  The  optimal  control 
u(t)  therefore  has  the  constant  gain  configuration  indicated  in  Figure  2. 
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More  generally,  if  equations  (6)  have  no  solution,  J(u(kT))  nay  be 
minimized  directly.  A  necessary  and  sufficient  condition  for  its  uncon¬ 
strained  minimization,  since  it  is  convex  in  u(kT),  is  7J(u(kT))  •  0.  This 
condition  reduces  simply  to 

(C0  ♦  D)TQ(C0  ♦  D)u(kT)  »  (CQ  ♦  D)TQ[z((k*l)T)  -  C4x(kT) }  ($) 

where  Q  -  (q_)  with  q.^  *  (6^  *  1,  i  *  j,  6..  *  0  otherwise) . Solutions 

of  (9)  can  be  written  in  the  form  (S).  If  equations  (6)  are  solvable,  then, 
in  view  of  the  convexity  of  J(u(kT))»  the  solution  spaces  of  (6)  and  (9) 
coincide.  Obviously  then,  attention  may  be  restricted  exclusively  to 
equations  (9). 

When  unobservable  states  are  present,  an  asymptotically  stable  digital 
observer  x  defined  by 

x((k*l)T)  x(kT)  ♦  Qu(kT)  ♦  G[y(kT)  -Cx(kT)],  k=0,l,2,...,  (10) 

is  employed.  In  this  case,  y  may  be  controlled  most  directly  by  minimizing 
J(u(kT))  under  the  assumption  of  complete  observability,  and  then  substituting 
the  estimate  x(kT)  for  x(kTl  in  (8).  Alternatively,  the  observer  output 
?((k+l)T)  •  Cfc((k+1)T)  may  be  substituted  for  y((k+l)T)  in  (5),  with 
x((k+l)T)  given  by  (10).  The  altered  J(u(kT))  may  then  be  minimized  as 
before.  The  optimal  control  again  assumes  the  basic  form  (8),  but  with  x(kT) 
replacing  x(kT)  and  an  additional  term  Hy(kT)  present  due  to  the  dependence 
of  y((k»l)T)  upon  y(kT).  Obviously,  the  effectiveness  of  either  of  these 
control  schemes  depends  upon  how  rapidly  x(kT)  converges  to  x(kT),  as  k 
increases,  for  given  initial  conditions  x(0)  and  x(0). 

Remarks 

3. I  Exact  output  matching  can  in  general  be  achieved  only  when  l  <  n,  that 


is,  only  when  there  are  at  least  as  many  control  variable  degrees  of  freedom 

as  there  are  controlled  variables.  In  practice,  if  Jt  »  m,  one  expects 

(6)  to  possess  a  unique  solution,  implying  that  C0  ♦  D  is  invertible,  and 

hence  that  c(kT)  =  0  in  (8).  If  l  <  m,  one  expects  that  the  values  of  m  -  i 

control  variables  may  be  assigned  freely,  so  that,  in  general,  c(kT)  +  0. 

Similarly,  multiple  solutions  in  (9)  are  expected  when  8,  <  m,  and  a  unique 

T 

solution  when  l  >  m.  The  later  implies  that  the  matrix  (CG  +  D)  Q(CG  +  D) 
is  invertible,  clearly  a  desirable  property. 

3.2  Because  J(u(kT))  is  continuous,  it  will  always  possess  a  constrained 
minimum  if  constraints  ai  <  u^(kT)  <  bi  are  imposed.  This  minimum  can  be 
determined  numerically  on  each  sampling  interval  using  standard  nonlinear 
programming  algorithms.  As  in  regulator  problems,  constraints  may  also  be 

~  T 

imposed  through  the  minimization  of  J(u(kT))  =  J(u(kT))  +u(kT)  R  u(kT),  R  >  0. 
Equating  Vj(u(kT))  to  zero,  a  system  of  linear  equations  for  the  components 
of  u(kT)  results  with  solutions  again  of  form  (8).  Alternatively,  constraints 
may  be  effectively  imposed  in  some  cases,  as  when  transferring  outputs  to 
set  points,  by  selecting  only  those  paths  z(t)  which  require  bounded  control 
energy  for  close  tracking. 

3.3  The  derivations  of  equations  (6)  and  (9)  are  based  simply  upon  the 

difference  equations  (3)  which  solutions  of  (1)  with  controls  (2)  satisfy. 

These  equations  are  valid  more  generally  if  A  »  A(t)  and  B  *  B(t)  are 

matrices  of  measurable  functions  which  are  bounded  on  finite  intervals  (see 

12 

Coddington  and  Levinson  ).  In  this  case,  however,  ♦  •  $(k)  and  6  *  0(k) 
depend  upon  k.  Consequently,  the  coefficient  matrices  in  (6)  and  (9)  and 
hence  the  gains  F  ■  F(k)  and  G  *  G(k)  in  (8)  also  depend  upon  k.  The  on-line 
implementation  of  the  control  scheme  discussed  here  may  then  present  serious 
computational  difficulties.  This  is  particularly  true  if  T  is  very  small. 


since  the  matrices  of  integrals  4(k)  and  9(k)  must  be  computed  and  a  system  of 
linear  equations  (6)  or  (9)  must  be  solved  at  each  sampling  instant  t  *  kT. 
Moreover,  if  unobservable  states  are  present,  the  additional  problem  of 
choosing  gains  G  *  G(k)  so  that  x  in  (10)  is  asymptotically  stable  must  be 
addressed. 


4.  NUMERICAL  EXAMPLES. 

Example  1.  A  model  foliating  problem. 

Let  a  one  dimensional  plant  be  described  by  the  state  and  output 
equations 


x(t)  ■  -x(t)  +  u(t) 
y(t)  *  x(t)  . 


(11) 


Let  the  equations  of  a  second  reference  system  be 

q(t)  *  [  2  l]  qCt)  *  [2]  r(t>  '  (12> 


Consider  the  problem  of  finding  a  control  u(t)  of  form  (2)  that  forces  y(t) 

to  match  z(t)  *  q^Ct)  closely  at  sailing  instants. 

»T 

In  (6),  C9  ♦  D  «  I  e"(^“T)<*T  ■  i  -  e“^  so  that  exact  output  matching 
'0 

can  be  achieved  using  u(t)  determined  by 

(13) 


u(kT)  «  Pz((k*l)T)  -  Gx(kT) 


9 

Figure  3  shows  z(t)  and  the  piecewise  exponential  function  y(t)  as  con¬ 
trolled  by  (13)  and  (14).  Here  T  =  1,  r(t)  is  a  unit  step  function,  and 
x(0)  a  q(0)  =  0. 

If  only  q2(t)  were  observable  in  (12),  it  would  be  necessary  to  equip 
q(t)  with  an  observer 

q((k*l)T)  =  0qq(kT)  ♦  0qr(kT)  +  Gq[q2(kT)  -q2(kT)]  , 

where 

.371074  .4444761  f.  6289261  _  f.  289579 " 

*q  =  ,  0q  =  ,  and  Gq  = 

888951  -.073402  _J  |_ .  8889S2  _[  [_.297672_ 

(See  Kuo1  for  digital  observer  design.)  Letting  z(t)  »  q^(t),  y(t)  could 
then  be  matched  to  qj(t)  using  (13)  and  (14).  A  control  simulation  in  this 
situation  is  pictured  in  Figure  4.  T,  r(t),  and  x(0)  and  q(0)  are  as  above. 
The  simulation  begins  with  q(0)  =0.  At  t  s  3  it  is  assumed  that  the  actual 
state  q(t),  evolving  according  to  (12),  is  subjected  to  a  disturbance  which 
.^sults  in  qj(3)  =  qj(3~)  ♦  .2  and  q2(3)  *  q2(3~)  +  -3.  Figure  4  illustrates 
the  convergence  of  y(t)  to  qj(t)  as  q(kT)  converges  to  q(kT)  from  q(3)  =q(3  ). 

Since  J(u(kT))  is  convex  in  a  single  variable,  it  is  easily  minimized 
subject  to  the  constraint  a  <  u(kT)  <  b.  Clearly,  if  VJ(a)  *  J'(ct)  *  0, 
to  minimize  J(u(kT),  choose  u(kT)  *  a  if  a  <  a  <  b,  and  choose  u(kT)  *  a 
if  a  <  a  and  u(kT)  >  b  if  a  >  b,  Figure  5  shows  the  example  of  Figure  4 
simulated  subject  to  the  constraint  0  <  u(kT)  <  1. 

Example  2.  A  digital  redesign  problem. 

o 

Kuo,  Singh,  and  Yackel  consider  the  problem  of  converting  a  continuous- 
data  control  system  into  a  sampled- data  system  while  preserving  the 
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performance  characteristics  of  the  original  system.  Sample  and  hold  devices 
are  inserted  into  the  continuous  system  and  its  input  and  feedback  gains  are 
modified.  The  modifications  insure  that  linear  combinations  of  the  plant 
states  of  the  resulting  sampled-data  system  match  the  same  combinations  of 
states  of  the  original  system  exactly  at  all  sampling  instants. 

The  output  matching  technique  of  section  3  may  also  be  applied  to  the 
digital  redesign  problem.  To  illustrate  this,  consider  the  simplified 
one-axis  Skylab  satellite  control  system  of  reference  9  with  plant  equations 

u(t)  ,  q (0)  =  0  (15) 

and  control 


0  1  “ 

'  o 

q(t)  = 

q(t)  + 

l 

0  0 

970,741 

u (t)  =  E(0)  r(t)  -  G(0)  x(t).  (16) 

V- 

The  scalar  input  r(t)  is  taken  to  be  a  unit  step  function,  and  the  input 
and  feedback  gains  E(0)  and  G(0)  are 

E(0)  =  11,800  and  G(0)  =  [11,800  151,800]  . 

Given  a  1  x  2  matrix  C,  formulas  for  modified  gains  E(T)  and  G(T)  are  derived 
in  reference  9  with  the  property  that  if  q^(t)  is  the  solution  of  (15)  with 
u(t)  given  by 

u(t)  =  E(T)  r(kT)  -  G(T)  x(kT),  kT  <  t  <  (k+l)T,  k  =  0,1,2,..., 

then  Cq^(t)  and  Cq(t)  match  exactly  at  all  sampling  instants. 

Let  x(t)  represent  the  solution  of  (15)  if  u(t)  is  of  form  (2),  while 


q(t)  is  the  solution  with  control  (16).  The  difference  equations  (3) 
satisfied  by  x(t)  are 
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Given  a  1  x  2  matrix  C,  in  the  framework  of  section  3,  the  redesign  goal  is 
to  determine  a  control  (2)  which  insures  the  close  matching  of  y(t)  =  Cx(t) 
and  z(t)  =  Cq(t).  By  (6),  if  the  scalar  C0  f  0,  y(t)  and  z(t)  may  be  matched 
exactly  at  sampling  instants  using  (13)  with  gains 

F  =  (CQ)_1  and  G  =  (C0)_1C$  .  (18) 

This  is  the  same  result  derived  in  reference  9  for  the  matching  of  Cq^(t) 
and  Cq (t) . 

Reference  9  discusses  in  detail  the  cases  C  =  [1  0]  and  C  =  [0  1], 
with  T  =  2.  The  corresponding  gains  (18)  for  these  cases  are 


C  -  [1  0] 
C  -  [0  1] 


and  G 


^  k 1  c  ‘ [1 01 

[o  ]  c  =  [0  1]  , 


(19) 


T 

where  0  ■  [0^  ©2]  •  Since  both  redesign  strategies  provide  exact  output 
matching,  the  simulation  results  for  the  above  choices  of  C  are  identical 
for  each  method.  Figures  4-13,  reference  9  (or  Figures  6.19  -  6.24,  Kuo*) 
illustrate  these  results.  Although  both  states  of  (15)  are  matched  well 
for  each  choice  of  C,  it  is  seen  that  better  simultaneous  matching  of 
states  occurs  when  C  ■  [01]. 

Kuo,  Singh,  and  Yackel  consider  the  digitalization  of  the  control  of 
a  continuous  plant  (1)  in  the  case  £  *  m.  (C  here  is  the  matrix  H  of 
reference  9.)  More  generally,  the  output  matching  technique  presented  here 
applies  to  the  digital  redesign  problem  when  £  and  m  are  arbitrary.  As 
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indicated  in  section  3,  when  D  =  0,  a  practical  condition  permitting  exact 
output  matching  is  the  existence  of  (C0)'1.  This  is  the  same  condition 
given  in  reference  9  for  the  case  l  ®  m. 

To  illustrate  the  more  general  problem,  letting  C  *  I,  the  identity 
matrix,  consider  the  Skylab  satellite  with  z(t)  *  Cq(t)  •  q(t).  For 
control,  minimize  the  weighted  sum 

J(u(kT))  *  jqJqjCCk+iyr)  -  x^k+UT)]2  ♦  £ q2[q2((k+l)T)  -x2((bl)T)]2. 
From  (9),  gains  F  and  G  for  the  optimal  control  (13)  are  given  by 

F  =  a_1(C0)TQ  and  G  =  o_1(C0)TQC*  (20) 

—  T 

where  Q  *  (q^)  has  9^  *6^  an<*  t^ie  scaler  a  a  (C0)  Q(C0)  is  assumed 
nonzero.  Figures  6,  7,  and  8  illustrate  simulations  of  q(t)  and  x(t)  as 
controlled  through  (13)  and  (20)  with  T  *  2,  =  2,  and  q2  *  1.  As  can  be 

seen,  close  matching  of  both  states  simultaneously  is  provided.  (The 
results  are  effectively  unchanged  when  *  1  and  q2  *  2,  or  when  q  ^  *  q  2  *  1 . ) 
In  fact,  when  the  actual  data  is  analyzed,  slightly  better  simultaneous 
matching  is  provided  by  C  !  I  than  by  either  C  =  (1  0]  or  C  *  [0  1]. 

To  illustrate  the  incorporation  of  the  digital  observer  of  Figure  2 
into  the  control  scheme,  consider  the  problem  of  matching  q2(t)  and  x2(t) 
when  only  x^(t)  is  observable.  Equip  x(t)  with  an  observer  (10)  with  gain 
G  *  [2  1/T]T.  For  control,  use  (13)  and  (19)  with  C  *  [0  1],  but  with  x(kT) 

replacing  x(kT)  in  (13).  A  control  simulation  under  these  conditions  is 
illustrated  in  Figures  9,  10,  and  11.  Here  again  T  *  2  and  q(0)  «  x(0)  »  0. 
It  is  assumed  that  5^(0)  *  0,  but  that  the  initial  estimate  x2(0)  *  .1  is  in 
error.  Since  x(kT)  effectively  converges  to  x(kT)  after  2  sampling  periods, 
x(t)  and  q(t)  are  closely  matched  for  t  >  6. 


mmmm 
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Example  3.  A  flight  control  example. 

In  Reference  4,  Mehra,  (et  al)  applies  the  Model  Algorithmic  Control 
(MAC)  technique  to  the  attitude  control  of  a  hypothetical  missile.  The 
three  axis  model,  with  independent  pitch  axis  and  coupled  roll-yaw  dynamics, 
is  described  by  the  state  equations 
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The  output,  state,  and  control  variables  are 


(21) 


y}  *  x,  *  angle  of  attack 

x2  =  perturbed  pitch  rate 
y2  *  x^  =  sideslip  angle 

*  perturbed  roll  rate 
Xg  *  perturbed  yaw  rate 
y^  =  =  roll  angle 

Angles  are  measured  in  degrees. 

With  the  missile  flying  at  Mach 


Uj  =  elevator  angle 
u2  =  aileron  angle 
Ug  =  rudder  angle 


at  20,000  ft.  and  weighing  239.5  lb.. 


the  parameters  in  (21)  have  values 
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It 

-1.4868 

Ng  =  290.48 

v  * 

2174.21 

Ma 

-149.93 

0^  =  .15708 

*1  = 

0.0 

V 

-.91237 

L,  =  8770.6 
<5p 

%  “ 

-281.11 

S  3 

-1559.2 

g  =  32.0 

N6r 

281.11 

Under  these  conditions,  consider  the  problem  of  transferring  the  outputs 
y-!  (t ) »  y2(t)  and  y3(t)  smoothly  from  given  initial  values  to  constant  set 
points  Cj,  c2«  and  c^  along  desired  trajectories  z^t) ,  z2(t),  and  Zj(t) 
respectively.  Letting  C  be  defined  so  that  y(t)  =  Cx(t),  the  matrix  C©  is 
invertible.  By  (6)  then,  exact  matching  of  y(kT)  and  z(kT),  k  =  1,2,..., 
may  be  achieved  using  u(t)  defined  through  (13)  with 

“-.846826  0.0  0.0 

F  =  0.0  .067512  5.269492  (22) 

_0.0  -.938841  . 080139  _ 

and 

“-.223562  -.060463  0.0  0.0  0.0  0.0 

G  -  0.0  0.0  -.107461  .001085  .003071  4.769538  (23) 

_0.0  0.0  .014207  .000023  .052151  .080906 

As  a  particular  example,  let  z^(t)  be  defined  at  sampling  instants  by 

^(kT)  *  ay.((k-l)T)  ♦  (l-o)c.,  i  -  1,2,3,  k»  1,2,... ,  (24) 

for  constant  a,  0  <  a  <  1.  (The  output  matching  technique  requires  only 
that  z(t)  be  specified  at  t  *  kT.)  Figures  12  and  13  show  a  simulation  of  (21) 

as  controlled  through  (13),  (22),  and  (23)  with  desired  trajectories  (24)  and 
set  points  Cj  ■  15",  c2  ■  10*,  and  ■  0*.  Here  x(0)  «  0,  T  »  .1,  a  »  .5, 
and  control  action  begins  at  t  ■  .5  (i.e.,  z(kT)  ■  0,  and  hence  u(kT)  »  0, 


Is 

0  <  k  <  5,  and  z(kT)  is  given  by  (24)  for  k  >  6). 

The  rapid  oscillations  of  u(t)  and  y(t)  every  .1  seconds  apparent  in 
Figures  12  and  13  are  clearly  undesirable.  Smoother  control  of  (21)  can 
be  realized  if  store  than  one  future  value  of  y(t)  is  considered  in  each 
control  computation.  As  an  example,  suppose  that  at  t  *  kT,  u(t)  is  to  be 
held  constant  at  value  u(kT)  for  tuo  sampling  periods.  Consider  the  mini¬ 
mization  of  the  two  step  ahead  matching  criterion 
3  3 

J(u(kT))«if  [zj((k*l)T)  -y .((k*l)T))2*i  l  fz?((k*l)T) -y .((k+2)T)]2  (25) 

Zi«l  X  1  ^i»l  X  1 

where  for  j  *  1,2,  z?((k*l)T)  *  z^((k+l)T)  as  defined  through  (24)  with  a  »  , 

Oj  >  a2*  Since 

y.((k*j)T)  =  CeA(jT)x(kT)  +cf  eA^T'TiBu(kT) ,  j  -  1,2, 

• 0 

J(u(kT))  may  be  minimized  by  solving  the  linear  equations  VJ(u(kT))  *  0. 

The  solution  of  these  equations  may  be  placed  in  the  form 

u(kT)  »  Flzl((k+l)T)  ♦  F2z2((k*l)T) -Gx(kT)  (26) 

1  2 

for  appropriate  gains  F  ,  F  ,  and  G.  Simulations  of  the  control  of  (21) 

using  controls  u(t)  determined  sequentially  through  (2)  by  the  u(kT)  in  (26) 

(i.e.,  u(kT)  is  applied  on  only  one  sampling- interval)  reveal  a  marked 

smoothing  of  u(t)  and  y(t).  ((25)  effectively  imposes  rate  contraints  upon 

the  components  of  y(t)  as  well  as  positional  constraints.)  The  simulation 

with  ot^  *  .75,  a2  »  .5,  and  x(0),  T,  and  the  c^  as  previously  specified  is 

pictured  in  Figures  14  and  15.  As  can  be  seen,  the  components  of  u(t)  and 

y(t)  converge  quickly  and  smoothly  to  steady  state  values.  The  gains 
1  2 

F  ,  F  ,  and  G  in  (26)  for  this  example  are 
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"-.110352 

0.0 

0.0 

"-.285081 

0.0 

0.0 

-1 

F  a 

0.0 

.127214 

1 . 135150 

F2  * 

0.0 

-.175194 

2.164551 

_0.0 

-.239291 

-. 002336 _ 

J 

% 

_0.0 

-.408877 

. 004776 _ 

and, 


"-.104394 

-.028234 

0.0  0.0 

0.0 

0.0 

G  * 

0.0 

o.o 

-.022853  .000639 

-.000910 

2.799690 

(28) 

0.0 

0.0 

.363652  .00001S 

. 0075S7 

.002551 

• 

R. 

K.  Mehra 

and  others 

(see  references  4,  5, 

11,  and  13) 

have  employed 

output  matching  extensively  in  the  control  of  linear  plants  with  impulse 
response  representation 

l  N 

y.  ((k+l)T)  *  1  l  a(n)..u.((k*l-n)T),  i  =  1,2,.  ...I  .  (29) 

1  j*l  n*l  3 

Letting  z-^kT),  k  *  1,2,...,  be  a  discrete  exponential  path  to  the  set 
point  Cj,  the  basic  impulse  response  technique  matches  yi ( (k-*-l)T)  to 
zi((k*l)T)  using 

u(kT)  *  K  lv(k) ,  v  (k)  -  z ( (k+1 )T)  -  £  \  a(n) . .u. ((k+l-nJT)  (30) 

j*l  n*2  11  3 

where  K  ■  (k^)  ^as  *jj  *  a(l).j.  Figures  33  a,  b,  reference  4,  show  a 
simulation  of  the  control  of  (21)  using  this  method  with  z^(kT) 

«  e‘V((k-l)T)  ♦  (1-e’1)^,  Cj  *  10°,  c 2  *  15*,  Cj  *  0*.  The  state  space 
model  (21)  was  used  to  generate  the  impulse  responses  a(n)y 

Figures  33  a,b  and  Figures  12  and  13  here  demonstrate  essentially 
equivalent  control  effectiveness  since  the  desired  output  trajectories  are 
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similar  for  both  examples  and  each  method  provides  exact  output  matching. 

(The  rapid  output  oscillations  of  Figures  12  and  13  are  absent  in  Figures 
33  a,b.  Evidently,  linear  interpolation  has  been  used  to  approximate  y(t) 
between  the  data  points  generated  by  (29)  and  (30).  Figures  12  and  13  were 
constructed  using  20  data  points  in  each  sampling  interval,  thus  revealing 
the  intersample  behavior  of  the  continuous  plant  (21).)  From  a  computational 
point  of  view,  however,  the  state  space  approach  offers  the  following 
advantages.  Each  control  computation  in  (13)  with  gains  (23)  requires  two 
matrix-vector  multiplications  and  a  vector  subtraction  using  the  stored 
gains  (system  constants)  F  and  G.  Each  computation  in  (30)  requires  a  series 
of  scalar  multiplications  and  additions  (to  determine  v(k)),  followed  by  a 
matrix-vector  multiplication,  using  the  system  constants  a(n)jj  and  K-1  and 
past  control  actions  u^ ((k+l-n)T) ,  i,j  *  1,2,3,  n  =  2,3,...,N.  Table  1 
compares  the  computational  requirements  of  the  two  control  methods. 

Methods  System  Constants  Algebraic  Operations/u(kT) 

State  Space  m(£  +  n)  *  27(15)  m(2£  ♦  2n  -  1)  *  51  (27) 

Impulse  Response  (N-l)i2  ♦  m2  *  450(250)  t2(2N-3)  +  i  ♦  m(2fc-l)  =  891(495) 

Table  1  -  Missile  Control  Computational  Requirements. 

(Since  the  value  of  N  employed  in  reference  4  is  not  specified,  it  is  assumed 

that  N  *  50  as  in  reference  11.)  Due  to  the  decoupling  of  states  x^  and  x^ 

1  2 

from  Xj  -  x6  in  (21),  the  matrices  F  and  G  in  (22)  and  (23)  and  F  ,  F  ,  and  G 

in  (27)  and  (28),  and  correspondingly  the  matrix  K*1  in  (26),  are  naturally 

partitioned.  Taking  this  into  consideration, the  actual  computational 
requirements,  given  in  parentheses  in  Table  l,  are  the  same  as  the  total 
requirements  for  two  and  four-dimensional  problems  considered  separately. 
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(E.g. ,  3(3+6)  ■  27  while  1(1+2)  ♦  2(2+4)  *  15.)  As  can  be  seen,  the  state 
space  approach  requires  significantly  fewer  parameters  to  define  control  with 
correspondingly  fewer  computations  per  control  calculation.  (Given  arbitrary 
plant  equations  (1)  with  A  ■  a  and  corresponding  impulse  response  representa¬ 
tion  (29) ,  the  formulas  in  Figure  16  remain  in  general  valid.  It  is 

2  2 

interesting  to  note  that  m  (2m  +  2n  -  1)  <  m  (2N  -  3)  ♦  m  ♦  m(2m-l)  whenever 
®  >  IjjTj  •  For  N  =  50,  this  implies  that  the  state  space  technique  requires 
fewer  operations  per  control  computation  roughly  as  long  as  m  >  ,  that  is, 

as  long  as  there  is  at  least  one  output  for  every  50  states.) 

As  a  final  comparison,  in  order  to  produce  control  smoothing,  reference 
4  introduces  matching  at  multiple  sampling  instants  through  control  blocking. 

A  gradient  search  algorithm  is  applied  to  the  nonlinear  optimization  problem 
which  results.  Simulation  results  using  this  approach  are  shown  in  Figures 
34,  37,  and  38,  reference  4.  These  figures  show  smoothing  comparable  to  that 
of  Figures  14  and  15.  The  minimization  of  (25)  through  (26),  however, 
yields  improved  system  performance  without  introducing  additional  computa¬ 
tional  complexity. 

V.  CONCLUSIONS 

A  simple,  direct,  and  quite  general  output  matching  approach  to  multi- 
variable  linear  digital  control  has  been  presented.  Digital  inputs  which 
force  plant  outputs  to  closely  track  desired  reference  trajectories  are 
generated  sequentially  by  controls  with  constant  forward  and  feedback  gain 
configuration.  Optimal  gains  are  determined  using  only  elementary  results 
from  linear  algebra  and  linear  systems  theory.  As  a  result,  the  control 
technique  discussed  is  characterized  by  its  conceptual  and  practical  simplicity. 
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Applications  of  output  matching  to  problems  in  model  following,  digital 
redesign,  and  direct  digital  design  have  demonstrated  its  flexibility,  ease 
of  application,  and  effectiveness.  Control  effectiveness,  of  course,  depends 
upon  sample  rate  and  desired  output  trajectory  selection.  An  intersample 
smoothing  strategy  promises  to  provide  improved  performance  without  increased 
sampling  rate  or  additional  computational  expense.  Since  each  linear  plant 
possesses  its  own  specific  dynamical  characteristics,  it  is  felt  that  in 
application,  only  experimentation  will  ultimately  reveal  those  best  choices 
of  sampling  rate  and  reference  trajectories  for  given  standards  of  acceptable 
system  performance. 
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Figure  3.  Y(t)  and  q,(t).  Example  1. 


Figure  7-  q _(t)  and  x,(t),  Example  2;  average  matching  of  states,  q  *  2,  q,  -  1 


Figure  11.  Control  functions.  Example  2;  q2(t)  and  x2(t)  matched,  x2(t)  unobservable. 
Continuous  control  u^(t)  given  by  (16),  digital  control  ud(t)  defined  through  (13)  and  (19) 


